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Abstract 

We present a series of observables for soft inclusive physics, and utilize them for comparison 
between two recently developed colour reconnection models; the new colour reconnection model in 
Pythia and the DIPSY rope hadronization model. The observables are ratios of identified hadron 
yields as a function of the final-state activity, as measured by the charged multiplicity. Since both 
considered models have a nontrivial dependence on the final-state activity, the above observables 
serve as excellent probes to test the effect of these models. Both models show a clear baryon 
enhancement with increasing multiplicity, while only the DIPSY rope model leads to a strangeness 
enhancement. Flow-like patterns, previously found to be connected to colour reconnection models, 
are investigated for the new models. Only Pythia shows a pj_-dependent enhancement of the A/K 
ratio as the final-state activity increases, with the enhancement being largest in the mid-p_L region. 
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I. INTRODUCTION 


The first run of the LHC has provided a large number of measurements probing both soft 
and hard QCD, and thereby a large number of tests for the Monte Carlo event generators. 
Even though the overall performance of the event generators have been quite good, there 
are still some phenomena that are insufficiently understood [1]. One of the more intriguing 
soft QCD deviations is the observed enhancement of A production [2, 3]. No model has 
been simultaneously able to describe the identihed hadron spectra at both LEP and LHC. 
This has led to the development of several phenomenological models [4-6], partly aimed to 
address this problem. With the planned low pile-up runs at the beginning of the second 
LHC run, it is now an ideal time to test these models further, and thereby probe the physical 
origin of the A enhancement. In this study we consider two of the models: the new colour 
reconnection (CR) model in the Pythia event generator [5, 7] and the colour rope model in 
the DIPSY event generator [4, 8, 9]. The models have previously been compared to pp data 
at i/i of 200, 900 and 7000 GeV. In this paper new possible observables to test the models 
are suggested, and predictions are made for collisions at i/i = 13 TeV. The observables are 
not model dependent, and can be used for constraining predictions from other models of 
soft inclusive physics. Both considered colour reconnection models are built upon the Lund 
model for string hadronization [10]. Nonperturbative differences can therefore be ascribed 
to differences in the new phenomenological ideas. 

One of the key ideas for the two models in question is jet universality. Stated in terms of 
the string model, it essentially means that fragmentation of a string does not depend on how 
the string is formed. Free strings at both lepton and hadron colliders should thus hadronize 
in a similar fashion. Fragmentation parameters are therefore tuned in the clean e’''e“ —)■ 
Z —)■ gg environment, and then directly applied to hadron colliders. Any discrepancy has to 
be due to physical phenomena not active at lepton colliders. For all the models attempting 
to describe the A enhancement, the enhancement is linked to the increased density of quarks 
and gluons in the hnal state at hadron colliders[11]. It would therefore be of natural interest 
to measure the A enhancement as a function of this density. The quark-gluon density is 
experimentally ill-dehned, however, and we suggest to use the number of charged tracks in 
the forward region as a measure of hnal-state activity. A similar idea for using the hyperon- 
to-meson ratio to search for indications of a miniQGP was suggested in ref. [12]. We suggest 
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ratios that allows for separation of strangeness enhancement from baryon enhancement, 
which both could be present in the hyperon-to-meson ratio. 

Another puzzling observation is the indication of collective effects in high-multiplicity pp 
collisions [13, 14], often interpreted as the presence of flow. These effects were only expected 
in the dense medium of heavy ion collisions, where the pressure gradients give rise to flow 
effects. A study of the models for pp collisions showed that CR generated similar effects 
even without the introduction of a thermalized medium [15]. We therefore consider one of 
the standard observables in heavy ion physics, that of identihed particle ratios as a function 
of px, separated into bins of centrality, and compare the model predictions for pp collisions. 
Since centrality is not experimentally well dehned in pp collisions, the number of charged 
tracks in the forward region is used as a measure of activity. 

The outline of the paper is as follows. In section II we will briefly recap the most important 
features of the two models considered. Comparison to existing data at = 91.2 GeV, 
and pp data at y/s of 200 GeV and 7 TeV, is shown in section III. The event selection and 
tuning for 13 TeV is described in section IV. In section V, the predictions at i/i = 13 TeV, 
for the second run of LHC, are presented. Finally, in section VI, we summarize and present 
an outlook. 


II. THE MODELS 

Both models for colour reconnection are built upon the Lund string model for hadroniza- 
tion. In this model, outgoing partons are connected with stringlike colour helds, which 
fragment into hadrons when moving apart. The model contains two main parameters rele¬ 
vant to this study, which determine the suppression of strange quarks and of diquarks (giving 
baryons) in the break ups. Assuming jet universality, these parameters are tuned to LEP 
data. 

Baryons can in addition be created around string junctions, which can arise as a con¬ 
sequence of colour reconnection. Consider the simple conhguration of two qq dipoles in 
hgure 1, which for example could have originated from a decay of two hF-bosons in a LEP 
environment, as described in ref. [16]. What essentially could be described as a quadrupole 
conhguration is instead described as either the original (on top) or the left conhguration 
in hgure 1. Without CR only the original conhguration is considered. Extending this type 
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FIG. 1. Sketch of how two qq dipoles (top) can be reconnected to different colour topologies (left 
and right). The right connection gives rise to a double junction, which in turn will produce baryons. 
Notice that the placement of the pairs differs in the junction figure. 

of colour reconnection to hadron colliders has been shown [17] to be a necessary condition 
to describe the rising of {pi_){Nch) distribntions. The QCD e-tensor gives rise to the right¬ 
most conhgnration, containing two jnnction connections, depicted as empty circles. Since 
snch junctions constitute proto-baryons, in the same way string segments constitute proto¬ 
mesons, they become an additional source of baryons. 


A. Colour reconnection in Pythia 

The new CR model in Pythia is situated just prior to the hadronization. It takes the 
leading-colour [Nc —)■ cxd) strings and transform them to a different colour configuration 
based on three principles: firstly the SU(3) colour rules from QCD determine if two strings 
are colour compatible {e.g. there is only a 1/9 probability that the top configuration of 
hgure 1 can transform to the left conhgnration purely from colour considerations). Secondly 
a simplistic space-time picture to check causal contact between the strings. Finally the A 
measure [18] (which is a string-length measure, A = Z]jlog(l + ^i/(2?Tio)) where the sum 
goes over all dipoles, rrii is the invariant mass of the dipole and tuq is a parameter) to decide 
whether a possible reconnection is actually favoured. Since the model relies purely on the 
outgoing partons, it is in principle applicable to any type of collision. So far it has only been 
tested for pp [5] and ee collisions [19]. The main extension compared to the other CR models 
in Pythia is the introduction of reconnections that form junction structures. From a pure 
colour consideration the probability to form a junction topology is three times larger than 
an ordinary reconnection. The junction will introduce additional strings, however, and it is 
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therefore often disfavoured due to a larger A measure. Given the close connection between 
junctions and baryons, the new model predicts a baryon enhancement. It was shown to 
be able to simultaneously describe the A production for both LEP and LHC experiments, 
which neither of the earlier Pythia tunes have been able to. 

The new CR model essentially contains two new parameters: a parameter that constrains 
the overall strength of the CR, and a parameter that controls the baryon enhancement. Both 
of these parameters were tuned to data [3, 20] from the LHC experiments at 7 TeV. 


B. Rope hadronization in DIPSY 

Rope hadronization [21] is a normer for QCD inspired models, which includes interactions 
between strings. From previous attempts to include this effect in Monte Carlo generators 
[22], it is well known that strange and baryonic content will rise in very dense events. 

A model introducing rope hadronization was recently developed and implemented in 
the event generator DIPSY [4]. Along with a hnal-state swing, the model introduces local 
calculation of string density, and corrects the evolution of the hnal-state parton shower and 
hadronization based on this local density. 

The model is based upon the idea [21] that when several parton pairs are next to each 
other in geometric space, they can act together coherently to form a colour rope. Each string 
is treated as a hux tube with a hxed radius, and the amount of overlap between strings, in 
impact parameter space and rapidity, can be directly calculated. 

If such an overlap is found to exist, it can have different effects, determined by SU(3) 
colour rules. The overlapping strings can end up in a colour singlet conhguration. This is 
handled by a hnal-state ’’swing”, that reconnects colour dipoles, in the hnal-state parton 
shower as the transformation from the top to the bottom left conhguration in hgure 1. In all 
other cases, the strings end up forming a ’’rope”. This is hadronized with a higher ehective 
string tension, rehecting the fact that more energy is available for the fragmentation, in 
accordance with results from lattice QCD [23]. In some cases, the strings forming the rope 
end up in a junction structure. In such cases the junction pair is handled using a simple 
approach, where the two junctions collapse to either two diquarks or two quarks, with a 
probability controlled by a tuneable parameter. The resulting strings are then hadronized 
with the appropriate ehective string tension. 
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An increased string tension results in more strange quarks and diquarks produced in 
string breakups. Since the effect increases with the density of quarks and gluons in the hnal 
state, the expected outcome is more baryons and strangeness among the resulting hadrons. 
The model includes two free parameters; the string radius and the probability for a junction 
to resolve to diquarks. Both are tuned to LHC data [3] at 7 TeV. 


III. COMPARISON TO DATA 

The models performs as intended when comparing to existing data. Ratios of baryons to 
mesons are enhanced for both models, whereas ratios of particles with strange quark content 
is enhanced only in the DIPSY rope model. Comparisons are done to ratios of integrated yields 
of identihed particles, using the analyses published through the Rivet [24] framework. The 
raw results from comparing the Monte Carlo to data using Rivet, are integrated to give 
hgures 2 and 3, using Matplotlib [25]. Error estimates are conservative, as they assume the 
error of all bins are fully correlated. 

In hgure 2, a comparison to LEP data [26-28] is seen. Two conclusions can be drawn 
from this figure. First of all, these are the data the original string model is tuned to. The 
fact that the Monte Carlo is so well aligned with data is thus not an indication that the 
string model predicts all these ratios so well, but rather that the parameters of the model 
are tuned to these data. The exception here is the H baryon[29], reflected in the H/S ratio, 
which lies below the observed value. However, the experimental statistical uncertainty is 
large for this ratio. 

The other conclusion, which is the most relevant for this article, is that only small effects 
at LEP data is observed. The A/A ratio increases by about 10 % for both the DIPSY rope 
model and the new CR model in Pythia (over their respective default models), but all 
models stay within the experimental uncertainty. The overall low variance is exactly what 
is expected, due to the low hnal state activity at LEP. 

In hgure 3 comparison to STAR data [30-32] at 200 GeV indicated that the description 
of the baryon to meson ratios improves with both models, while the description of the S/A 
ratio only improves with the DIPSY rope model. The change in the A^/tt ratio is not visible 
on this scale for this energy. 

Comparison to 7 TeV data from ALICE [33, 34] and CMS [3] conhrms that the description 
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FIG, 2, Comparison to e+e at 91,2 GeV from ALEPH, SLD and PDG, Color online. 
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FIG, 3, Comparison to pp data at 200 GeV from STAR (left) and at 7 TeV from ALICE and CMS 
(right). Color online. 


improves, even for the fl/S ratio. The description of the p/n ratio is seen to be somewhat 
worse with the new models. This could either have a mundane explanation originating in 
the fact that the very low-p_L area of the individual distributions (which contains most of 
the multiplicity) are not fully understood, or have further reaching consequences. We point 
to the measurements suggested in the next section of this paper to shed light on this issue. 


7 




































































































IV. TUNING AND EVENT SELECTION 


Before studying exclusive observables at 13 TeV, it is necessary to verify that the baselines 
for the two models agree reasonable well. Normally this is achieved by tuning the models 
to the available data. Data at = 13 TeV is, however, not yet published in a state 
where event generators can be tuned to it, so the DIPSY model was instead tuned to the 
Pythia predictions for diVc/i/d?7, {p±) {Nch) and the multiplicity distribution. Both models 
will eventually have to be retuned, when more data, in a suitable format for tuning, become 
available. Only small effects are expected from the retuning, firstly due to fragmentation 
mainly being determined from LEP data, and secondly since the already presented resnlts 
at 13 TeV show a good agreement between the Monash tnne and the data [35, 36]. The fnll 
list of all parameters changed from their defanlt valnes is inclnded in an appendix. 

An event and particle selection was implemented to mimic a possible experimental setup. 
Each particle is required to have p± > 0.15 GeV. Two different p regions are used; a 
forward region (2 < |? 7 | < 5) to measure the activity, and a central region (|? 7 | < 1) to 
measnre the identihed hadron yields. The reason for the split is to avoid any potential 
bias, which otherwise happens at low Nch, in particular for ratios involving both charged 
and non-charged hadrons. Since DIPSY does not have a model for diffraction, only non- 
diffractive events are considered for both models. To reflect this in the event selection, only 
events with at least six forward charged particles are considered. 

All particles with cr > 10 mm are treated as stable. In practice this means that 
71^, K, A,'E and D are all stable whereas 0 (which decays strongly) is not. This intro- 
dnces some donble connting in the 0/A-ratio, where a 0 can potentially be connted in the 
nnmerator and its decay prodncts in the denominator. 


V. PREDICTIONS FOR 13 TEV 

Differences between the colonr reconnection models are best determined using observ¬ 
ables controlled by hadronization effects. Ratios of identihed particles is exactly snch an 
observable, since particle species prodnction is determined by the qnark and diqnark content 
in string breaks. In the hrst part of this section, ratios of identihed particles are shown as a 
fnnction of Nch in the forward region, as a measnre of event activity. Then how-like ehects 
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FIG. 4. Ratios of identified hadrons as functions of at ^/s = 13 TeV. The top row shows 

meson ratios with the numerator having one more strange quark than the denominator. The middle 
row shows baryon to meson ratios, with same amount of strange quarks. The bottom row shows 
baryon ratios with the numerator having one more strange quark than the denominator. Note that 
the vertical axis differs between the figures and that zero is suppressed. 

are considered, by showing [A/K) {p±) in fonr different bins of Nch in the forward region. 


A. Particle ratios 

Ratios of hadrons with different strange and baryon nnmbers as fnnction of event activity, 
measnred as fnnctions of are shown in fignre 4. The strangeness enhancement in meson 

prodnction is probed by the K/n and ratios, for which the nnmerator always contains 
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one more strange quark than the denominator. As expected, only the DIPSY rope model 
shows an enhancement relative to the baseline, since it contains a strangeness enhancement. 
The new Pythia CR model lies slightly below the baseline. This can be explained by phase- 
space constraints for low invariant-mass strings, which the new model produces more of. It 
should be recalled that both the new as well as the old models are capable of describing the 
total Kg yield at 7 TeV. Thus, the limited effects in this ratio is somewhat expected. The 
(p/K ratio shows more promise as a means to distinguish between the two models, since the 
DIPSY model shows a larger enhancement. It is, however, more experimentally challenging. 

The baryon enhancement is tested for both hadrons containing zero or one strange quark, 
p/tt and KjK. For both ratios, and both models, clear enhancements are expected and seen. 
For the KjK ratio both models agree quite well, which is not surprising, given that both 
models are tuned to describe the inclusive KjK distributions at 7 TeV. A similar picture is 
seen for the p/tt ratio, indicating similar predictions for the baryon enhancement from both 
the models. 

The multistrange baryon enhancement is tested in the same way as the strange-meson 
enhancement by considering the ratios S/A and hl/S. The large variations at low multiplicity 
for both distributions are due to statistical fluctuations. For S/A the DIPSY rope model shows 
a clear enhancement as opposed to the new Pythia CR model. The A/p ratio is not shown, 
but the enhancement is similar to the enhancement of S/A. An enhancement is seen for both 
models in the hl/S, with the enhancement factor being around 2.5 for the DIPSY rope model 
in the highest multiplicity bins. This is larger than any of the other enhancements seen. 
The enhancement for the new Pythia CR model is somewhat surprising, as the increased 
junction production should be equal for both S and hi. The production of hi in the standard 
Pythia fragmentation is, however, signihcantly suppressed, as the production of ss-diquarks 
is disfavoured. This suppression is not present in the junction handling, since it takes two 
already formed quarks and combine into a diquark. The enhancement in the new Pythia 
model should therefore not be interpreted as a ’’real” strangeness enhancement, but more 
as an absence of suppression of ss diquarks. For the DIPSY model the above effect is also 
present, but there is an additional enhancement of strangeness and diquarks. It should be 
noted that the hi baseline from LEP is not that well constrained, due to a large experimental 
uncertainty, and the model predictions are below the actual measurements. A measurement 
of (fl/S) [Kch) would cast light on whether an actual activity-based enhancement takes 
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place. 


Increased hyperon production in high activity pp events have previously been associated 
with production of a miniQGP [12], The hyperon-to-pion ratio is only indirectly shown in 
hgure 4, but the rise is similar to the one predicted by miniQGP. The new models therefore 
provide an alternative explanation, if such an enhancement is observed. 


B. Flow-like effects 

The K/K ratio as a function of p± for different ranges is shown in hgure 5. The 

two models show different behaviours for the different multiplicity ranges: the DIPSY rope 
model only gives a small enhancement (~ 10% at maximum) between the lowest and highest 
multiplicity regions. Even though the differential enhancement is generally below 10 %, the 
enhancement of the ratio of integrated yields is about 20 %, which is in good agreement 
with hgure 4. It should be noted that the DIPSY model is inadequate in describing the high 
p_L tails (p_L >~ 4 GeV). This was observed for 900 GeV and 7 TeV in ref. [4]. 

The new Pythia GR model shows a clear change in px with increasing multiplicity. 
The enhancement is largest in the mid-p_L region {p± ~ 2 — 6 GeV), leading to a ’’peak” 
structure. This structure looks qualitatively similar to what is observed in PbPb and pPb 
collisions [37, 38]. The peak also moves towards larger p± with increased multiplicity, an 
ehect normally attributed to radial how in heavy ion collisions [39]. That the new GR model 
predicts a qualitatively similar ehect in pp collisions is quite intriguing and strengthens the 
hint at a potential connection between how and GR ehects already observed [15]. 


VI. CONCLUSIONS 

A series of new model-independent observables, well suited for distinguishing between 
between diherent physical models for soft inclusive physics is suggested. The observables 
are ratios of identihed hadrons measured as a function of event activity, with the identihed 
hadrons chosen such that a distinction is made between baryon-only, strangeness-only and 
baryon-and-strangeness enhancement. Measurement of these observables at present and 
future energies at pp colliders is encouraged, as the can serve as constraints on any soft 
physics model aiming to explain low-multiplicity and minimum bias data simultaneously. 
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FIG. 5. Ratio of A/X as a function of p± in three bins of In the right column the new 

colour reconnection models are shown, and in the left column the old ones. 

The observables are, in this article, used to separate two new CR models. The new CR 
model in Pythia only contains a baryon enhancement with increasing multiplicity, while 
the DIPSY rope models contains both a baryon and a strangeness enhancement. The multi- 
strange hyperon ratios, as well as the (f)/K ratio, provide clear observables for distinguishing 
between the two models. It should be mentioned that this is already possible to observe in 
inclusive measurements, but the separation into different multiplicity regions highlights the 
enhancement. 

Both new models are based on interactions between strings in the hadronization phase, 
and confirmation of the common predictions made by the two models is a direct hint that 
colour reconnections among strings are of physical importance. Both baseline models show 
almost no dependency on multiplicity for the identified hadron yield ratios. Therefore, any 
observed dependency would provide a clearer indication that the old models miss a feature. 
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better than an inclusive measurement alone could provide. We therefore strongly suggest 
that these observables should be measured at the LHC experiments. In this paper we only 
studied the effects at a center-of-mass energy of 13 TeV, but the effects should also be visible 
in the already collected data at 7 TeV. 

We have also shown that one of the CR models predicts effects similar to those normally 
attributed to radial flow in heavy ion collisions. This is in agreement with earlier indications 
that also hint at a connection between the two phenomena. It should however be recalled 
that neither of the models provide a satisfactory description of the individual p± spectra for 
the identihed hadrons. And before these are fully understood, claims of connections between 
flow and CR may be premature. 
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Appendix A: Model parameters 

A complete list of all the parameters that differ from their default values for the considered 
models. 

Hadronization model parameters are found in table I, Pythia parameters in table II and 
DIPSY parameters in tables III and IV. 
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Pragmentation parameter 

Pythia def. Pythia new 

DIPSY DIPSY rope 

StringPT:sigma 

0.335 

0.335 

0.32 

0.31 

StringZ:aLund 

0.68 

0.36 

0.30 

0.41 

StringZ:bLund 

0.98 

0.56 

0.36 

0.37 

StringFlav:probQQtoQ 

0.081 

0.078 

0.082 

0.073 

StringFlav:ProbStoUD 

0.217 

0.22 

0.22 

0.21 


0.5 

0.0275 

- 

- 


0.7 

0.0275 

- 

- 

StringFlav:probQQltoQQOj oin 

0.9 

0.0275 

- 

- 


1.0 

0.0275 

- 

- 


TABLE I. Table of parameters of the string hadronization model, which differs from Monash tune 
[35] default values. The changed parameters have been retuned to LEP and SLD data, cf. figure 2 


Pythia parameter Default New 


MultiPartonInteractions:pT0Ref 2.28 2.15 


BeamRemnants:remnantMode 0 1 
BeamRemnants:saturation - 5 
ColourReconnection:mode 0 1 
ColourReconnection:allowDoubleJunRem on off 
ColourReconnection:mO - 0.3 
ColourReconnection:allowJunctions - on 
ColourReconnection:junctionCorrection - 1.2 
ColourReconnection:timeDilationMode - 2 
ColourReconnection:timeDilationPar - 0.18 


TABLE 11. The new Pythia CR model introduces a number of new parameters, and requires 
retuning of a few old ones, besides hadronization. The details of the retuning can be found in ref. 


[5]. 
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DIPSY parameter 

Default Rope 

FragmentationScheme default dipole 

StringRO 

0.773 

StringmO 

0.113 

BetaPopcorn 

0.2 


TABLE III. The DIPSY rope model introduces three extra parameters, which are fixed using pp 
data from LHC. See ref. [4] for the meaning of the parameters. 


Energy pp 200 GeV pp 7 TeV pp 13 TeV 


DIPSY parameter Default Rope Default Rope Default Rope 


LambdaQCD 

0, 

.29 

0, 

.26 

0, 

.17 

0. 

.25 

0, 

.29 

0. 

.27 

RMax 

2, 

.32 

3, 

.34 

3, 

.23 

2. 

.90 

1, 

.05 

3. 

.39 

PMinusOrdering 

1, 

.05 

0, 

.98 

1, 

.24 

0. 

.67 

0, 

.31 

0. 

.75 

PTScale 

0, 

.70 

0, 

.92 

1, 

.60 

1. 

.65 

1, 

.28 

1. 

.35 


TABLE IV. The DIPSY initial state model needs retuning at each energy to reproduce total charged 
multiplicity. See ref. [9] for the meaning of the parameters. 
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